Novel Quantum Spin Assisted Tunneling in Half Metallic Manganite Tunnel Junctions 
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The electron tunneling in half metallic manganite tunnel junctions is studied by using a quantum 
mechanically treated double exchange model. We show that the stimulation of spin excitations, 
caused by the strong Hund's coupling between the conduction eg electrons and the localized quantum 
spins of the manganite ions, would assist electrons to tunnel through the junction even with the 
antiparallel aligned magnetizations of the electrodes. This mechanism gives rise to an extra tunnel 
conductance, in addition to that predicted by Julliere's model. Our theory is in good agreement 
with the voltage dependence of the tunnel conductance in manganite tunnel junctions observed by 
experiments. 

PACS numbers: 71.30.+h, 72.15.Rn, 71.28.+d, 71.27.+a 



Tunnel magnetoresistance (TMR) in ferromagnetic 
junctions (FJ), first identified more than two decades 
ago 10], is of both fundermental interest and potential 
applications to magnetic sensors and memory devices. A 
simple model Q proposed by JuUiere is widely adopted 
to account for the observed TMR, in which the conduc- 
tance of the polarized electrons are proportional to the 
product of the density of states (DOS) at the Fermi level 
from both magnetic electrodes for each spin channel. Ac- 
cording to this model, the TMR ratio is given by 
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where Rp and Rap are the tunnel resistances when the 
magnetizations of the two ferromagnetic electrodes are in 
parallel (P) and antiparallel (AP) alignments. P here is 
the spin polarization of the electrodes. This model works 
quite well for Fe, Co and Ne based ferromagnets. From 
Eq. (|I]) the TMR ratio tends to become infinity when the 
electrodes are made of half-metallic ferromagnets where 
only one spin component exists and P = \, such as in 
some mixed valence manganites Lai^ajAajMnOa, (A=Ca, 
Ba, Sr, etc.). The bulk of these manganese oxides, well 
known as the colossal magnetoresistance (CMR) materi- 
als, have been intensively studied in recent years. The 
involved physics is based on the double exchange (DE) 
model which can be described by the following 

Hamiltonian 



(2) 



where the two terms represent the conduction Cg elec- 
tron hopping between neighboring Mn ions and the on- 
site Hund's coupling between the electron spin a and the 
localized t2g core spin Si of Mn'^^ . It is well known that 
in mixed valence manganite system the strong Hund's 
coupling Jh » t dominates the basic physics by sep- 
arating the two spin subbands. The even stronger on- 
site Coulomb repulsion (not written out in Eq. (^)) fur- 
ther enhances the energy difference between the two sub- 
bands. As a result, a reasonable assumption of a full spin 



polarization (P = 1) was extensively accepted in the lit- 
eratures. For this reason, the manganese oxides have 
been recognized as being good candidates for the study 
of various spin polarized tunnclings |^-^. The authors 
in Refs. |^-|| fabricated FJ using Lao.erSro.asMnOs, 
Lao.8Sro.2Mn03 and Lao.rSro.sMnOa to obtain a high 
TMR ratio 80%, 150%, and 430% respectively [|-|. 
More recently, a TMR ratio as high as 730% are ob- 
tained by Jo et al. with Lao.rCao.sMnOs as electrodes 
1^. All these results were explained by Julliere's model, 
from which a spin polarization P = 0.54,0.65,0.83 and 
0.88 are respectively assumed. While these values are 
much larger than that of a typical common magnetic 
transition metal, they are still far away from the expected 
P = 1 for half-metallic manganites. Furthermore, it is 
found that the tunnel conductance G in these mangan- 
ite junctions depends on the bias voltage V. In both 
the P and AP magnetization alignments, G has a con- 
siderable increase when V rises up to tens of mV 
As a result, TMR is strongly reduced. Bratkovsky at- 
tributed the y-dependence of G to the lowering of the 
barrier height of the insulator layer by the applied volt- 
age, but his calculation shows that this effect becomes 
appreciable only when \eV\ > 0.3eV jl^. Therefore the 
origin of such a voltage dependence of the conductance 
can not be explained by the existing theories ||l|,|l0| and 
so far still is a puzzle. To account for the discrepancy be- 
tween the expected full spin polarization of the observed 
finite TMR ratio, it was proposed that the presence of 
defect states in the barrier or resonant state can reduce 
the TMR 0. A reduced TMR was also obtained by 
Lyu et al. who introduced spin-flip matrix elements in 
the tunneling Hamiltonian . Recently, Itoh et al. also 
studied the effect of the temperature dependence of the 
spin polarization on TMR Jl^ . 

In this work, we investigate the voltage-dependent elec- 
tronic tunneling in the half-metallic manganite junctions 
based on a quantum mechanical treatment of the DE 
model. It is found that in the AP aligned electrodes 
even without any spin-flip scattering from the tunnel- 
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ing Hamiltonian, an electron can still tunnel through the 
junction resulting in a nonzero conductance and thus 
finite TMR. Such a tunneling comes from the intrinsic 
quantum spin effect, and is accompanied by the stimula- 
tion of spin excitations in the ferromagnetic electrodes. 
In the classical spin 6*^00 limit, our result reduces to 
a formula similar to that of JuUiere's model. For finite 
S, the quantum spin effect becomes prominent and sup- 
presses strongly the TMR amplitude. We show that our 
theory explains well the voltage dependence of the tunnel 
conductance in manganite tunnel junctions observed in 
Refs. [| and g. 

We start from the following model for manganite tun- 
nel junctions 



H = Hr 



(3) 



where and Hji are the quantum DE Hamiltonian with 
the form of Eq. (||) in the left and right electrodes. Ht 
represents the spin-conserved incoherent tunneling pro- 
cess of conduction electrons through the insulator barrier. 
In the momentum space, it is given by 



(4) 



According to the perturbation theory, the net elec- 
tronic tunnel current is 



m,n pp' ,cr 

X (e-^E^ - e-^''-)S{Er, - + eV) , (5) 



where |m), \n) are eigenstates of Hl + -ff_R with energy 
and En, Z =Tre-''(^^+^«), with (3 = {keT)-^ the 
inverse temperature, and V is the applied bias voltage. 
The delta function ensures the conservation of energy in 
the tunneling process. 

To derive the tunnel current from Eq. (||), it is con- 
venient to use the Schwinger-boson or slave-fermion rep- 
resentation 14 1 in which the electron operator Cia is 
expressed as a combination of a spinless charge fermion 
fi and a neutral spin boson bi„. In the present strong 
Hund's coupling {Jh — s- 00) case, the spin of an electron 
on a site is aligned parallelly to the local spin, and the 
electron operator can be written as — ff"bf^ / ^/2Si + 1 



for i £ L, and 



1,/VWTT for 



i <E R, where Si — 3/2 is the local spin at the site i, 
and [/(JO-' (^) are matrix elements of the 2x2 matrix 
U{e) = cos(6'/2)i -t- isin{e/2)ay. Without loss of gen- 
erality, we assume the magnetization in the left electrode 
is parallel to the z-axis while the magnetization in the 
right electrode makes an arbitrary angle 9 with the z-axis. 
9 — Q and tt correspond to the P and AP configurations 
of the junction magnetizations. The total spin of a site 
(the local spin plus the spin of the itinerant electron) 
is 5*1°* = \^„„ib\c,(ya<y'biai , with the local constraint 



J2a^iabi<j = 25'; + f} fi- We use a mean-field approxi- 
mation that was derived in Ref. for the unperturbed 
Hamiltonian {a — L, R) 



E 
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where ojk ~ pk^, with the spin stiffness p — po{l — 
Efe7fe[exp(/3po7fe) - ll^VlS-S}, and = ^(7*. - 6) with 
t = tj:^{blM/2S. Here po = <(.///,)/25, 7fe - 
6 — 2(cosA:2: -t- cosky + coskz) and S = Si + {I — x)/2 
is the average spin per site. The total number of spin 
bosons is n-f -|- — 2NS, where n„ = J2k(^la^ka) 
and N is the total number of sites. At temperatures 
below the Curie temperature Tc and at each electrode, 
there are macroscopic numbers of bosons condense into 
k = and a =t state, no.| = — = 2N Sm, where 
m ~ (rij — +ni) = M/Ms, with M the magne- 

tization of the electrode and Mg its saturation value, is 
just the normalized magnetization of the electrode. Em- 
ploying the Schwinger boson representation for Eq. (j^) 
and after some algebra, we obtain the following tunnel 
current 
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where nl 



-<^f) + and ' 



(e 



-q — \.'^ < ■''ko 

the distribution functions of charge fermions and spin 
bosons. In the derivation we have treated |Tpp'P as a 
constant \T\'^. Taking into account the condensation of 
the spin bosons in the fc = and cr —'\ state, and that the 
band width of the fermion is in the order of electron volts 
much larger than energy scales of ksT and eV, thus the 
DOS of fermions near the Fermi level is approximated as 
a constant, the dynamical conductance of the junction 
G = dl/dV can be derived as 
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G (9) = |Go{l + COS0 + + '-^ J gb{oj)doj 



e\V\ 



1 



J+e\V\ 



gb{uj)dujj^ [n''iu;)-n''{^')]gbico')dLo'} , (8) 

where Go - ^{2§rfi)^\^\^9j{eF) with gf{e) the DOS 
of fermions, n''{Lj) = {e^"^ — oJm — {6n'^)^^^p is the 

maximum energy of the spin excitation, and gbi^^) is the 
normalized DOS of bosons, i.e., " gb{(^)duj = 1. 

Eq. (H) is the tunnel conductance of the half-metallic 
manganite for arbitrary magnetization configuration, in 
which a voltage dependence is found in the last two 
terms. All the last three terms are related to the spin 
size 5" and vanish when the spin is treated classically 
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{S oo). This indicates that they come from the quan- 
tum spin effect. In the classical spin limit the TMR re- 
duces to AR/R = 2m^/{l — rn^), which has a similar 
form to JuUiere's formula Eq. (|l]), but the spin polar- 
ization P there is replaced by the normalized magnetiza- 
tion m here. P and m are usually different, the former 
is the relative difference of the DOS at the Fermi level 
between the two electron subbands split by the Hund's 
coupling, in the large Jh limit it is always 1 (see, e.g., 
the calculation of P in Ref. jl2j) so the system is half- 
metallic, while the latter reaches 1 only at zero temper- 
ature, and decreases with the increase temperature due 
to the spin wave excitations. When the quantum spin ef- 
fect is taken into account, even at zero temperature with 
full magnetization (m=l), the fourth term contributes a 
nonzero conductance. In the application of a small volt- 
age e\V\ < LJrn, this term reduces to (e|y|/p)^/^/67r^S', 
so that our theory predicts the tunnel conductance G{d) 
at low temperatures behaves as A + where A, B 

are constants independent of V . In Fig. 1(a) we show 
a comparison between our theory and experimental data 
in Ref. in the AP magnetization configuration of the 
Lao.eTSro.asMnOa/SrTiOa/Lao.eySro.asMnOa junction at 
T = 4.2K. The electronic hopping t is taken to be 0.35 
eV so that p ~ 14meV. Apart from the difference of 
a small voltage-independent constant, which may come 
from spin-flip centers, defect states in the barrier, m < 1 
at finite temperature or imperfect alignment of moments 
in a less-than-ideal sample |l^Jll[], our calculated con- 
ductance is in excellent agreement with the experimental 
measurement. For comparison the conductance in the 
classical spin treatment {S = oo) of the DE model is also 
plotted by the dashed line, which is independent of V. 
The dependence of the *^on(|^^|a^^c^ c)i^ S^ru^S indi- 
cates that we have to specify the value of the bias voltage 
when discussing the TMR. Interestingly, for the case of 
e\V\ larger than the band width of spin excitations, both 
the fourth and the last terms in Eq. (||) become inde- 
pendent of the bias voltage. While further increasing of 
the conductance with the rise of the voltage may be still 
possible through lowering the energy barrier of the in- 
sulator layer [1^, such an increase does not change the 
ratio of TMR. From Eq. (||) and under this large volt- 
age, the TMR ratio in the full magnetization m = 1 of 
an ideal half metallic manganite junction is AR/R = 2S, 
which is less than 400% for realistic S < 2. A divergent 
TMR ratio can be obtained only when S is treated as 
classically {S — > oo) so that the quantum spin effect in- 
duced tunneling vanishes and the conductance in the AP 
magnetization configuration goes to zero. 

The effect of temperature comes into the conductance 
mainly through the normalized magnetization to, which 
decreases with increasing the temperature. In Fig. 1(b) 
we plot the voltage-dependent conductance at differ- 
ent temperatures. To compare our calculations with 
the experimental results in Ref. [|| in which the FJ 



was Lao.8Sro.2Mn03/SrTi03/Lao.8Sro.2Mn03, we need 
to determine the magnetizations at different tempera- 
tures. Here the magnetizations are estimated from the 
measurement in Ref. p^ , they are taken as m = 1 at 
T = 0, TO = 0.6 at T 0.4Tc, and to = 0.2 at T = O.ST^, 
with Tc = 290K g. In Fig. 1(b) at each temperature 
we shift the voltage-independent constant conductance 
so that the comparison of the ^-dependent conductance 
at different temperatures can be shown more clearly. Fig. 
1(b) shows that the increase of the conductance with 
the increasing voltage becomes slower at higher temper- 
atures, in agree with the experimental results of Ref. 
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FIG. 1. Tunnel conductance G(7r) as a function of the 
applied voltage (a) at temperature T = 4.2K and doping 
X = 0.33, where the conductance has been shifted a small 
voltage-independent constant so as to compare with the ex- 
perimental measurement of Ref. Q, which is shown by the 
dots in the present figure, and Go in Eq. (^) is taken to be 
6.11 the dashed line is for S = oo, which is voltage 

independent, (b) at temperatures T = 0, 0.4Tc and O.STc and 
X = 0.2, with Tc = 290K. The electron hopping is taken to be 
t = 0.35eV. 

At this stage, we wish to address on why it is possible 
for electrons to tunnel through the junction for a spin 
conserved tunneling Hamiltonian when the magnetiza- 
tions of the two fully polarized electrodes are antiparallel 
to each other. Although the Schwinger boson method 
is convenient for calculation, the spin boson operator in 
it describes the total spin on a site as a whole, which 
conceals the actual motion of different spins. To see the 
physical origin of the nonvanishing G(7r), we now start di- 
rectly from the DE Hamiltonian Eq. (||) . At each site the 
Hund's coupling —J^Si ■ oi has two eigenvalues —JhSi 
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and Jh{Si + 1), where the former is 25"; + 2 fold degener- 
acy with total spin Si + 1/2, and the latter is 25"; degen- 
eracy with total spin 5; — 1/2. In the case of the strong 
Hund's coupling, at an occupied site only the lower en- 
ergy state can exist, i.e., the spin of the electron must 
be parallel to the local spin. In this case, the effective 
DE hopping of the eg electrons between two neighboring 
sites is usually treated to be te// = t cos 9ij where 
9ij is the relative angle between local spins on the two 
sites. When the effective hopping is generalized to man- 
ganite junctions with AP configuration at zero temper- 
ature, 9ij = IT for i and j in different electrodes, and 
the electron can not tunnel through the junction. That 
is exactly what JuUiere's model predicts. However, such 
a classical picture neglects the quantum nature of local 
spins. In fact, when the quantum effect is taken into ac- 
count, the spin of a conduction electron in one electrode 
can not be perfectly antiparallel to a localized spin in the 
other electrode. For example, when a spin-up {sz — 1/2) 
electron in the left electrode moves to a spin-down site 
{Sz = —Si) in the right electrode, the formed |l/2, —Si) 
state, which is not an eigenstate of the Hund's coupling, 
has nonzero components in both the spin parallel (the Cg 
electron spin and the local spin are in parallel) and spin 
antiparallel subspaces. Using the method of projection 
operator, it is straightforward to show 

where hAp)_ = (11/2,-5,) + V^i\ - 1/2, -Si + 
l))/^/257TTand|Vap) = iV^i\l/2,-Si)-\-l/2,-Si+ 
l))/y/2Si + 1 are states in which the electron spin and 
the local spin are in parallel and antiparallel. Both [ipp) 
and \ipap) are eigenstates of S^^^ (z-component of the to- 
tal spin of the site) with eigenvalues —Si + 1/2. At the 
same time, they are also eigenstates of the Hund's cou- 
pling — JnSi • (fi, with eigenvalues —JhSi for \'ipp) and 
JniSi + 1) for \-ipap)- For large Jh, state IV'ap) is for- 
bidden due to its very high energy. However, Eq. (|^) 
shows that |l/2,— 5'/) has a nonzero component in \ipp), 
i.e., there is a finite probability 1/(25/ -I- 1) for the spin 
of a spin-up electron to be parallel to a spin-down local 
site, so that even in this Jh — *■ oo case the electron may 
still tunnel through the junction. This tunneling process 
vanishes when Si oo, which is an implication that this 
effect is induced by the quantum nature of the localized 
spins. In the second term of \ipp), the spin of the electron 
is flipped while the local spin of the site is twisted. Such a 
spin-flip and spin-twist process is an intrinsic dynamic ef- 
fect originated from the Hund's coupling and the discrete 
quantum spin orientations. Moreover, the twist of the lo- 
cal spin can spread to the whole system to further lower 
the energy, thus stimulates a spin excitation and enables 
electrons to tunnel through the barrier. This mechanism 
is very different from that considered in some other pa- 



pers ||Tl|;|l^,|T^ . The energy scale involved in this tun- 
neling process is the spin excitation energy which comes 
from the applied voltage and thus causes the induced 
tunneling to be voltage dependent. 

Finally, it is worth to mention that in some experi- 
ments, at temperatures 77K and above, a quadratic vari- 
ation of the conductance with the applied voltage were 
found This is likely due to the fact that the DOS 

of the spin excitations near the ferromagnet/insulator 
(F/I) regions is different from that in the bulk man- 
ganites. Experimental study of the magnetic proper- 
ties at surface boundary of the half-metallic ferromagnet 
Lao.ySro.aMnOa indeed indicates that except at very low 
temperature (T < 30K), the temperature dependence of 
the surface-boundary magnetism is significantly different 
from that of the bulk ||l^ . In the present theory, a V"^ de- 
pendence of the conductance can be produced if the spin 
excitation spectrum has the form uj^ = + PJ-I^J-li 
where and k± are the wave vectors in and perpen- 
dicular to the plane of the F/I interface, respectively. 
Such a dispersion relation was proposed for the lightly 
doped Lai_a:Sra:Mn03 with layered antiferromagnetic 
spin alignment jl^ , and it may exist in the F/I inter- 
face regions of the half metallic manganite junctions. 
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